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SUMMARY 


OBJECTIVE 

Investigate  the  validity  of  a  parabolic  equation  (PE)  model  for  predicting  radio  field  strengths  in 
horizontally  inhomogeneous  environments  by  performing  comparisons  between  the  model  and 
experimental  data. 


RESULTS 

Excellent  agreements  were  found  at  VHP  and  UHF  frequencies  with  good  agreement  in  S-  and 
X-bands.  In  some  cases,  the  predicted  curves  for  the  S-band  comparisons  underestimated  that  of  the 
measured  data  at  large  ranges.  This  may  be  the  result  of  phenomena  such  as  surface  roughness, 
backscatter,  etc.,  not  accounted  for  in  the  model.  Discrepancies  may  also  result  from  the  presence  of 
evaporation  ducts  not  included  in  the  environmental  inputs  to  the  model  because  of  a  lack  of  detailed 
measurements.  This  would  account  for  lower  predicted  signal  levels  at  higher  frequencies. 


Accession  For  | 

NT  IS 

GRAS;  I 

OTIC  TAB 

□ 

U:)r’:-no';n'_'cd 

□ 

J','". ;  ' 

n  ... 

Bv 

L  i  :  t  r 

*  ^  0 ' 

/ 

A'/a' 

laSUl'. 

7  r.i'lea 

Av  jJ.  end /or  | 

Dlst 

ial 

th 

CONTENTS 


1.0  INTRODUCTION  .  1 

2.0  BACKGROUND  .  2 

2.1  Refractive  Etfects  .  2 

2.2  Split-Step  PE  model  .  6 

2.3  Propagation  Factor  and  Path  Loss  .  10 

2.4  Numerical  Implementation  .  11 

3.0  VALIDATION  OF  PE  MODEL  .  13 

3.1  Comparison  with  Range-Independent  Environments  .  13 

3.2  Comparison  with  Range-Dependent  Environments  .  14 

4.0  CONCLUSIONS  .  61 

REFERENCES  .  62 

RELATED  READING  .  63 

APPENDIX:  GUADALUPE  ISLAND  PROFILE  TABLES . A-1 

FIGURES. 

1.  Refractivity  N  and  modified  refractivity  M  versus  altitude  for  a  standard 

atmosphere .  3 

2.  Refractivity  N  and  modified  refractivity  M  versus  altitude  for  a  surface-based  duct 

created  by  an  elevated  trapping  layer  .  3 

3.  Raytrace  diagram  for  a  transmitter  at  25  meters  (m)  in  a  100-m  surface-based  duct 

resulting  from  an  elevated  layer .  4 

4.  Raytrace  diagram  for  a  transmitter  at  25  m,  standard  atmosphere  conditions  .  5 

5.  M-profile  for  a  typical  evaporation  duct .  5 

6.  Height-gain  plot  comparison  of  PEPC  and  MLAYER  for  300-m  surface-based  duct  ...  15 

7.  Height-gain  plot  comparison  of  PEPC  and  MLAYER  for  20-m  evaporation  duct  .  16 

8.  Measured  M-unit-versus-height  profiles  from  Canterbury  Project  environmental 

measurements  for  5  Augu.st  1947  .  16 

a 


9.  X-band  comparison  between  PEPC  and  Canterbury  Project  measurements  for 

5  August  1947  17 

10.  S-band  comparison  between  PEPC  and  Canterbury  Project  measurements  for 

5  August  1947  17 

1 1 .  Range-dependent  ray  trace  with  trapping  layer  (represented  by  shaded  area)  varying 

with  range.  Antenna  height  is  25  m  .  18 

12.  Coverage  diagram  for  environment  of  figure  11  .  19 

13.  Environment  and  radio  data  record  from  Guadalupe  Island  measurements  for 

12  March  1948  .  20 

14.  Environment  and  radio  data  record  from  Guadalupe  Island  measurements  for 

8  April  1948  21 

15.  Environment  and  radio  data  record  from  Guadalupe  Island  measurements  for 

13  November  1947  .  22 

16.  Interpolated  M-unit-versus-height  profiles  for  fixed  ranges  from  12  March  1948,  with 

variation  in  trapping  layer  height  and  thickness  represented  by  shaded  areas  .  23 

17.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  12  March  1948, 

along  slant  path  from  69  to  89  nmi .  23 

18.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  12  March  1948, 

along  slant  path  from  91  to  101  nmi  .  24 

19.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  12  March  1948, 

along  slant  path  from  1 15  to  130  nmi  .  24 

20.  Comparisons  between  PEPC  and  mea.sured  radio  data  at  520  MHz  for  12  March  194i<, 

along  slant  path  fri)m  69  to  89  nmi .  25 

21.  Compari.sons  between  PEPC  and  measured  radio  data  at  520  MHz  for  12  March  1948, 

along  slant  path  from  91  to  101  nmi  .  25 

22.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz  tor  12  March  1948, 

along  slant  path  from  1 15  to  130  nmi  .  26 

23.  Comparisons  between  PEPC  and  measured  radio  data  at  ?3U0  MHz  for  12  March 

1948,  along  slant  path  from  80  to  90  nmi  .  26 


Hi 


24.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  12  March 

1948,  along  slant  path  from  91  to  96  nmi  .  27 

25.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  12  March 

1948,  along  slant  path  from  124  to  130  nmi  .  27 

26.  Coverage  diagram  for  environment  of  12  March  1948  at  170  MHz .  28 

27.  Coverage  diagram  for  environment  of  12  March  1948  at  520  MHz .  29 

28.  Coverage  diagram  for  environment  of  12  March  1948  at  3300  MHz  .  30 

29.  Coverage  diagram  using  homogeneous  profile  measured  at  39  nmi  taken  from  the 

12  March  1948  data  record .  31 

30.  Comparisons  between  PEPC  and  measured  radio  data,  using  the  average  profile  from 

the  12  March  1948  data  record  for  170  MHz  .  32 

31.  Comparisons  between  PEPC  and  measured  radio  data,  using  the  average  profile  from 

the  12  March  1948  data  record  for  3300  MHz .  32 

32.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  8  April  1948, 

along  slant  path  from  60  to  72  nmi .  33 

33.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  8  April  1948, 

along  slant  path  from  81  to  95  nmi .  33 

34.  Comparisons  between  PEPC  and  mea.sured  radio  data  at  520  MHz  for  8  April  1948, 

along  slant  path  from  60  to  72  nmi .  34 

35.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  8  April  1948, 

along  slant  path  from  1 18  to  136  nmi  .  34 

36.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  8  April  1948, 

along  slant  path  from  97  to  107  nmi  .  35 

37.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz  for  8  April  1948, 

along  slant  path  from  81  to  95  nmi .  35 

38.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz  for  8  April  1948, 

along  slant  path  from  97  to  107  nmi  .  36 


39.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz  for  8  April  1948, 

along  slant  path  from  1 18  to  136  nmi  .  36 

40.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  8  April  1948, 

along  slant  path  from  60  to  69  nmi .  37 

41.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  8  April  1948, 

along  slant  path  from  85  to  95  nmi .  37 

42.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  8  April  1948, 

along  slant  path  from  97  to  104  nmi  .  38 

43.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  8  April  1948, 

along  slant  path  from  126  to  136  nmi  .  38 

44.  Coverage  diagram  for  environment  of  8  April  1948  at  170  MHz  .  39 

45.  Coverage  diagram  for  environment  of  8  April  1948  at  520  MHz  .  40 

46.  Coverage  diagram  for  environment  of  8  April  1948  at  3300  MHz  .  41 

47.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  13  November  1947, 

along  slant  path  from  41 .5  to  52.5  nmi  .  42 

48.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  13  November  1947, 

along  slant  path  from  56  to  75  nmi .  42 

49.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz  for  13  November  1947, 

along  slant  path  from  77  to  90  nmi .  43 

50.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz  for  13  November  1947, 

along  slant  path  from  41 .5  to  52.5  nmi  .  43 

51.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz  for  13  November  1947, 

along  slant  path  from  56  to  75  nmi .  44 

52.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz  for  13  November  1947, 

along  slant  path  from  77  to  90  nmi .  44 

53.  Ctimparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  13  November  1947, 

along  slant  path  from  41 .5  to  50  nmi  .  45 

54.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  13  November  1947, 

along  slant  path  from  61  to  75  nmi  .  45 

v 


55.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz  for  13  November  1947, 

along  slant  path  from  77  to  87  nmi .  46 

56.  Coverage  diagram  for  environment  of  13  November  1947  at  170  MHz  .  47 

57.  Coverage  diagram  for  environment  of  13  November  1947  at  520  MHz  .  48 

58.  Coverage  diagram  for  environment  of  13  November  1947  at  3300  MHz  .  49 

59.  Data  records  from  Canterbury  Project  for  19  June  1947  .  50 

60.  Data  records  from  Canterbury  Project  for  11  July  1947  .  51 

61.  Comparison  between  PEPC  and  measured  radio  data  at  3240  MHz  from  Canterbury 

Project  data  record  for  19  June  1947  .  52 

62.  Comparisons  between  PEPC  and  measured  radio  data  at  9875  MHz  from  Canterbury 

Project  data  record  for  19  June  1947,  with  antenna  height  at  26.5  feet  .  52 

63.  Comparisons  between  PEPC  and  measured  radio  data  at  9875  MHz  from  Canterbury 

Project  data  record  for  19  June  1947,  with  antenna  height  at  86  feet  .  53 

64.  Coverage  diagram  for  frequency  and  environment  from  figure  61  .  54 

65.  Cover  diagram  for  frequency  and  environment  from  figure  62  .  55 

66.  Cover  diagram  for  frequency  and  environment  from  figure  63  .  56 

67.  Comparisons  between  PEPC  and  measured  radio  data  for  flight  path  labeled  GH  from 

Canterbury  Project  data  record  of  1 1  July  1947,  for  antenna  heights  at  26.5  feet  .  57 

68.  Comparisons  between  PEPC  and  measured  radio  data  for  flight  path  labeled  GH  from 

Canterbury  Project  data  record  of  1 1  July  1947,  for  antenna  height  at  86  feet  .  57 

69.  Comparisons  between  PEPC  and  measured  radio  data  for  flight  path  labeled  HI  from 

Canterbury  Project  data  record  of  1 1  July  1947,  for  antenna  height  at  26.5  feet  .  58 

70.  Comparisons  between  PEPC  and  measured  radio  data  for  flight  path  labeled  HI  from 

Canterbury  Project  data  record  of  1 1  July  1947,  for  antenna  height  at  86  feet  .  58 

71.  Coverage  diagram  for  frequencies  and  environments  from  figures  67  70,  antenna 

height  at  26.5  feet  .  59 

72.  Coverage  diagram  for  frequencies  and  environments  from  figures  67-70,  antenna 

height  at  86  feet  . .  60 


V/ 


I.O  INTRODUCTION 


Many  field-strength  prediction  models  for  tropospheric  radiowave  propagation  over  the  ocean  have 
been  written,  based  on  the  limitation  that  the  environment  is  laterally  homogeneous.  It  was  found  that 
the  assumption  of  a  horizontally  stratified  troposphere  led  to  valid  propagation  assessments  86%  of  the 
time  [1],  While  the  troposphere  over  the  sea  does  exhibit  horizontal  homogeneity  over  relatively  long 
distances  in  most  cases,  meteorological  conditions  occur  occasionally  in  which  the  environment  may 
change  drastically  in  just  a  few  kilometers,  such  as  at  air-mass  boundaries  associated  with  wave  cyclones 
and  land/ocean  interfaces  [2],  It  is  safe  to  say  that  horizontally  varying  environments  do  occur  14%  of 
the  time.  This  paper  addresses  the  validity  of  the  parabolic  equation  (PE)  model  for  handling  these 
environments. 

The  more  conventional  method  of  modeling  electromagnetic  wave  propagation  in  the  troposphere 
i<i  performed  by  normal  mode  thet)ry,  but  most  modal  techniques  are  based  on  a  horizontally  stratified 
atmosphere.  A  parabolic  equation  method,  originally  developed  by  Fock  in  1946,  allowed  the 
environment  to  vary  with  range  as  well  as  height  [3]. 

Two  methods  can  be  used  to  solve  the  PE.  One  method  approximates  the  derivatives  in  the  PE  by 
finite  differences.  The  other  solves  the  PE  (with  some  approximations)  by  the  split-step  Fourier  method 
developed  by  Hardin  and  Tapper!  [4|.  The  split-step  method  requires  a  constant  vertical  mesh  size 
because  of  the  Fast  Fourier  Transforms  (FFT)  used,  which  may  require  large  data  arrays  for  sufficient 
sampling.  The  main  advantage  the  split-step  method  has  over  finite  difference  is  that  it  is  more  efficient 
at  higher  frequencies.  The  disadvantage  is  that,  with  higher  frequencies,  one  needs  to  have  a  larger 
transform  size  (smaller  vertical  mesh  size).  The  finite-difference  method  has  the  advantage  that  variable 
height  and  range  steps  can  be  used  to  keep  the  number  of  data  points  stored  at  a  minimum.  Over  the 
relatively  long  distances  at  which  one  is  interested  in  looking,  the  split-step  method  is  more  numerically 
efficient  because  of  the  FFTs  employed.  The  split-step  algorithm  thus  allows  real-time  predictions  based 
on  current  environment  descriptions. 

While  the  acoustic  community  has  been  using  PE  modeling  for  some  time,  it  was  only  recently 
(within  the  last  decade)  that  the  radar  community  has  applied  the  PE  method  to  tropospheric  radiowave 
propagation  [5-  8]. 

The  importance  of  the  split-step  PE  method  is  not  only  that  it  provides  an  exact  solution  to  field 
equations  (within  the  limits  of  the  operator  approximation)  for  a  homogeneous  atmosphere  in  a  relatively 
short  time,  but  that  it  can  also  predict  (with  relatively  small  errors)  field  strengths  for  laterally 
inhomogeneous  environments.  The  importance  of  taking  into  account  an  environment  that  varies  with 
..nge  will  become  apparent  as  case  studies  are  presented  between  experimental  data  and  the  predicted 
fields. 

A  parabolic  equation  model  for  the  personal  computer  (PEPC)  was  developed  by  Professor  Fred 
Tapper!  of  the  University  of  Miami  in  conjunction  with  the  Naval  Ocean  Systems  Center  (NOSC), 
San  Diego,  during  1989  and  1990.  Comparisons  between  predicted  PE  fields  and  measured  radio  and 
meteorological  profiles  on  over-water  paths  from  Guadalupe  Island  to  San  Diego  in  1947  and  1948  are 
presented  |9|.  Some  comparisons  also  were  made  again.st  experimental  data  taken  from  the  Canterbury 
Project  ( 10|.  The  appendix  contains  a  list  of  profiles  used  for  the  figures  presented  in  this  report. 
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2.0  BACKGROUIVD 


2.1  REFRACTIVE  EFFECTS 

The  refractive  index,  n,  of  air  is  defined  as  the  ratio  of  the  velocity  of  propagation  of  an 
electromagnetic  (EM)  wave  in  a  vacuum  to  that  in  air.  Since  EM  waves  travel  more  slowly  in  air  than 
in  a  vacuum,  this  ratio  is  always  greater  than  one.  At  the  Earth's  surface,  the  numeric  value  of  n  is 
usually  between  1 .0(X)250  and  1 .000400.  A  more  convenient  number  to  use  is  the  refractivity,  N,  which 
is  defined  as  N  =  (n  -  1)10^.  Normal  surface  N  values  then  range  from  250  to  400.  Refractivity  can 
be  expressed  as  a  function  of  atmospheric  pressure,  air  temperature,  and  humidity  as  follows; 

N  -  (77.6P)/T  +  (3.73*10^  e)/T- 

where  P  is  the  pressure  in  millibars,  T  is  the  temperature  in  degrees  Kelvin,  and  e  is  the  vapor  pressure 
in  millibars.  In  a  "standard"  well  mixed  atmosphere,  both  temperature  and  humidity  decrease  with 
increasing  height  such  that  N  decreases  at  the  rate  of  about  39  N-units  per  kilometer.  An  EM  wave 
propagating  initially  horizontal  to  the  Earth's  surface  will  be  refracted  downward,  but  with  a  rate  of 
curvature  less  than  that  of  the  Earth.  If  the  air  temperature  should  increase  with  height,  or  the  humidity 
decrease  abnormally  fast  with  altitude,  N  will  decrease  faster  than  normal  with  height.  If  N  decreases 
faster  than  157  N-units  per  kilometer,  an  EM  wave  will  be  refracted  downward  with  a  radius  of  curvature 
greater  than  the  Earth's  surface.  When  this  occurs,  the  EM  wave  can  be  trapped  in  a  surface  duct.  The 
EM  wave  is  first  refracted  downward  toward  the  surface  and,  if  the  surface  is  sufficiently  smooth,  the 
wave  is  specularly  reflected.  This  process  can  be  repeated  to  ranges  far  beyond  the  normal  horizon. 
Such  continuous  refractions/reflections  form  a  surface-based  duct  and  allows  radar  detections  at  greatly 
extended  ranges. 

A  more  convenient  method  for  expressing  the  refractivity  to  determine  the  occurrence  of  ducting 
is  with  the  use  of  modified  refractivity,  or  M-units.  M  is  defined  in  terms  of  N  as 

M  =  N  +  0.157  h 

for  h  in  meters.  Modified  refractivity  accounts  for  the  curvature  of  the  Earth  and  ducts  can  be  determined 
by  inspection  of  M  plotted  versus  height.  Whenever  M  decreases  with  height,  a  trapping  layer  is  formed 
and  indicated  by  a  negative  M-gradient.  Examples  of  N  and  M  plotted  versus  height  for  a  standard 
atmosphere  (118  M-units  per  kilometer)  are  given  in  figure  1.  In  this  figure,  M  increases  with  altitude 
and  no  trapping  layer  or  duct  is  formed.  Figure  2  shows  a  plot  where  M  decreases  with  height  and  forms 
a  trapping  layer  (in  this  case,  a  surface-based  duct).  If  the  M  value  at  the  top  of  the  trapping  layer  is  less 
than  that  at  the  surface,  then  a  surface-based  duct  is  formed.  If  the  M  value  at  the  top  of  the  trappin; 
layer  is  greater  than  that  at  the  surface,  an  elevated  duct  is  formed;  the  vertical  extent  of  the  duct  is 
determined  by  the  height  below  the  trapping  layer  where  the  M  value  is  equal  to  M  at  the  top  of  the 
trapping  layer.  Two  other  terms  are  used  to  describe  N-  or  M-gradients  other  than  standard  or  trapping, 
namely  subrefractive  and  superrefractive.  Superrefractive  means  more  bending  than  standard  refraction 
but  less  than  trapping.  Subrefractive  means  less  refraction  than  standard.  Those  atmospheres  generally 
do  not  have  the  spectacular  effect  on  EM  system  performance  that  ducting  does. 
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Figure  1 .  Refractivity  N  and  modified  refractivily  M  versus  altitude  for  a  standard  atmosphere. 


Figure  2.  Refractivity  N  and  modified  refractivity  M  versus  altitude  for 
a  surface-based  duct  created  by  an  elevated  trapping  layer. 
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Another  measure  of  refractivity  that  was  widely  used  in  the  1940s  is  the  B-unit,  defined  as 

B  =  N  +  0.039  h 

for  h  in  meters.  Thus,  for  a  standard  atmosphere,  B  would  be  invariant  with  height. 

The  two  distinct  types  of  surface  ducts  that  affect  naval  EM  systems  are  surface-based  ducts  from 
elevated  layers  and  evaporation  ducts.  Surface-based  ducts  from  elevated  layers  are  formed  by  a  sharp 
change  in  the  index  of  refraction  gradient  between  a  cool,  moist  marine  air  mass  and  a  dryer,  warmer  air 
mass  above  it.  These  ducts  generally  affect  all  systems  with  fre^iuencies  greater  than  about  100  MHz, 
provided  both  transmitter  and  receiver/target  are  in  or  near  the  duci.  These  ducts  tend  to  be  on  the  order 
of  100  to  300  meters  thick.  This  type  of  duct  is  quite  common  off  the  Soutiiern  California  coast,  where 
the  same  temperature  inversions  that  cause  smog  to  be  trappt^  in  the  Los  Angeles  air  basin  also  create 
strong  surface-based  ducts  at  sea.  The  duct  has  a  dramatic  effect  on  the  signal  level  well  '  'yond  the 
horizon  but,  generally,  does  not  greatly  affect  the  signal  level  within  and  near  the  horizon.  Quite  often, 
a  "skip  zone"  is  formed  where  ranges  just  beyond  the  horizon  are  in  the  normal  shadow  zone,  while 
slightly  greater  ranges  are  illuminated  by  energy  that  has  been  refracted  down  by  the  trapping  layer.  The 
range  to,  and  extent  of,  the  skip  zone  is  obviously  a  complex  function  of  the  duct  height,  the  M-unit 
gradients,  and  the  transmitter  height.  An  example  of  such  a  skip  zone  is  given  in  figure  3,  which  is  a 
raytrace  diagram  for  a  25-meter  transmitter  antenna  in  a  100-meter  surface-based  duct.  The  trapping  layer 
is  represented  by  the  shaded  area.  The  skip  zone  is  from  approximately  12  to  17  nmi.  For  comparison, 
the  raytrace  of  figure  4  is  at  the  same  transmitter  height  but  under  standard  atmospheric  conditions. 


Figure  3.  Raytrace  diagram  for  a  tran.smitter  at  25  meters  (m)  in  a 
100-m  surface-based  duct  resulting  from  an  elevated  layer. 
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Figure  4.  Raytrace  diagram  for  a  transmitter  at  25  m,  standard  atmosphere  conditions. 


Evaporation  ducts  are  created  by  a  rapid  decrease  in  humidity  from  the  air/sea  interface.  These 
ducts  are  limited  by  boundary-layer  mechanics  to  perhaps  50  meters,  although  heights  of  less  than  20 
meters  are  the  most  common.  Figure  5  is  an  example  of  a  typical  evaporation-duct  M-unit  profile.  The 
evaporation-duct  height  corresponds  to  the  minimum  on  the  M-unit  profile  and  is  a  measure  of  the 
strength  of  the  duct.  Unlike  the  surface-based  duct  resulting  from  an  elevated  layer,  the  radar  and  the 
target  do  not  have  to  be  "in"  the  duct  for  enhanced  signal  levels  to  be  observed.  These  ducts  primarily 
affect  EM  systems  with  frequencies  greater  than  3  GHz  and,  if  a  surface-based  duct  from  an  elevated 
layer  is  present,  it  will  be  the  dominant  propagation  mechanism  for  ranges  well  beyond  the  horizon. 


Figure  5.  M-profile  for  a  typical  evaporation  duct. 
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2.2  SPLIT-STEP  PE  MODEL 


A  brief  discussion  will  be  given  of  the  PE  method  and  how  the  program  calculates  field  strength. 
For  more  detailed  descriptions  of  the  theory,  see  (3,  4,  11],  In  the  following  discussion,  the  troposphere 
is  assumed  to  vary  in  range  and  height  only,  making  the  field  equations  independent  of  azimuth.  The 
Earth's  surface  will  be  considered  perfectly  conducting,  and  only  horizontal  polarization  will  be  addressed. 

As  with  every  electromagnetic  wave  problem,  one  begins  with  Maxwell's  equations  for  steady-state 
sinusoidal  (harmonic)  time-varying  fields.  Second-order  partial  differential  equations  (p.d.e.)  for  the 
electric  and  magnetic  fields  can  be  derived  from  the  coupled  first-order  differential  equations  by  the  usual 
method.  For  a  horizontally  polarized  electric  dipole,  the  only  nonzero  component  of  the  electric  field  is 
E(p  (r,0),  where  r,  0,  (p  represent  spherical  coordinates.  To  write  the  second-order  p.d.e.  in  the  form 
of  the  scalar  Helmholtz  wave  equation,  a  change  of  variable  is  used 

\lr(r,0)  =  7rsin0  E<p(r,0)  (1) 

and,  changing  from  spherical  to  rectangular  coordinates  by  using  z  =  r  -  a  and  x  =  a0  (applying  the 
Earth-flattening  transformation,  valid  for  z/a«:l,  x/a«l),  one  obtains  the  two-dimensional  elliptic  wave 
equation 


+ 


_  +  m(x,z) 


tlf(x,z)  =  0 


m(x,z)  =  n(x,z)  +  z/a 
kg  =  2nlX 


where  now 


(2) 


i|i(x,z)  =  \Jx  E(p(x,z)  (3) 

and  x  represents  the  horizontal  spatial  coordinate  (range),  z  is  the  vertical  coordinate  (height),  a  is  the 
Earth's  radius,  k^  is  the  vacuum  wavenumber,  and  m(x,z)  is  the  modified  index  of  refraction  taking  into 
account  the  Earth's  curvature  and  neglecting  small  terms,  since  we  are  only  interested  in  the  field  at  large 
ranges  (kxM). 

In  much  of  the  literature  that  has  been  published  on  the  theoretical  description  of  the  derivation  of 
the  PE,  the  parabolic  approximation  is  stated  in  two  ways.  The  first  method  removes  rapid  phase 
variations  (on  the  basis  that  one  is  only  interested  in  field  variations  that  are  large  compared  to  a 
wavelength)  by  performing  an  envelope  transformation 

i|/(x,z)  =  u(x,z)  e*  ”  (4) 


which  is  sub.stituted  in  equation  2  to  give 


az“ 


+  2i  k. 


dx 


2  7 

+  k„  (m^ 


1) 


u(x,z)  =  0. 


(5) 
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TTie  horizontal  variation  in  u(x,z)  is  now  assumed  to  be  slow 


a^u 

«2k„ 

au 

ax“ 

O 

ax 

(6) 


in  which  case  one  is  justified  in  neglecting  the  second-order  partial  derivative  in  the  horizontal  component 
in  equation  5.  The  second  method  is  to  define  an  operator  Q(x)  and  factor  equation  2: 


_a 

dx 


ikoQ(x) 


^  +  ikoO(x) 
ox 


Il/(x,z)  +  ikg 


4.QW 

ox 


iKx,z)  =  0 


(7) 


Q(x)  = 


1 


+  m‘ 


(8) 


where 


dx 


,Q(x) 


=  ^  -  0-. 
dx  dx 


(9) 


The  quantities  in  brackets  in  equation  7  represent  the  equations  for  outward  and  inward  propagating 
waves,  respectively.  Since  we  are  only  interested  in  outward  propagating  waves,  we  use 


dx 


ili(x,z)  =  ik^  Q(x)  1(1  (x.z). 


(10) 


After  substitution  of  equation  4,  the  final  equation  we  must  solve  is 


dx 


u(x,z)  =  ikj,  0(x)  u(x,z) 


(11) 


where 


0(x)  =  Q(x)-  1. 


(12) 
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This  is  now  reduced  to  an  initial  value  problem  that  neglects  backscattering  and  can  be  solved  by 
"marching  out"  in  range.  Equation  1 1  is  referred  to  in  the  literature  as  the  general  parabolic  equation 
(GPE),  and  is  exact  for  a  range-independent  environment. 

While  it  is  not  incorrect  to  make  the  assumption  given  by  equation  6,  the  assumption  necessarily 
restricts  accurate  solutions  to  the  field  to  near-horizontal  directions  or  small  propagating  angles.  The 
second  method  has  the  advantage  that  approximations  can  be  made  for  backscattered  waves;  the  different 
methods  of  approximating  or  "splitting"  the  operator  Q(x)  lead  to  different  parabolic  equation 
propagators. 

For  the  simplest  case,  assume  that  the  index  of  refraction  in  0(x)  is  range  independent  (m  is 
constant);  then  the  field  at  some  range  x  +  Ax  can  be  found  from  u(x); 

u(x  +  Ax)  =  u(x)  -  e‘'‘’'^^^’‘^(x).  (13) 


At  this  point,  Q(x)  can  be  "split"  to  get  ultimately  a  form  for  u(x  -l-  Ax)  suitable  for  numerical 
computation  on  a  computer.  Let  us  define  two  new  variables  a  and  b: 


1  .  2  , 
a  =  —  — ;  b  =  m  -  1. 
.2  ;^2 
ko  dz 


(14) 


Q(x)  can  now  be  written  as  a  function  of  a  and  b  as 


0(x)  =  (1  +  a  +  b)^^^  -  1. 


(15) 


The  two  most  widely  used  approximations  tor  Q(x)  are 


Q,(x)  =  (!/2)a  +  (l/2)b 

{16a) 

Q2(x)  =  (Ua)'^2  ^  „  _  1 

(16b) 

where  in  Qj  (x)  the  term  under  the  radical  in  equation  15  has  been  expanded  in  a  binominal  series  with 
higher-order  terms  dropped.  This  approximation  leads  to  the  standard  parabolic  equation  (SPE).  O2  (x) 
leads  to  the  wide-angle  parabolic  equation  (WAPE)  originally  derived  by  Feit  and  Fleck  [12]  to  study 
propagation  in  optical  fibers,  and  later  applied  to  acoustic  propagation  by  Thomson  and  Chapman  [13]. 

Thomson  and  Chapman  have  done  extensive  error  analysis  for  the  approximations  0i  (x)  and 
Q2  (x)  and  have  stated  that  Q2  (x)  gives  much  more  accurate  numerical  predictions  at  higher  angles. 
However,  their  application  was  toward  acoustic  propagation,  in  which  severe  ducting  conditions  (large 
gradients)  in  the  ocean  can  trap  modes  corresponding  to  relatively  large  propagation  angles.  For  the 
electromagnetic  case,  large  gradients  found  in  the  troposphere  are  two  orders  of  magnitude  smaller  than 
that  found  in  the  ocean.  So  while  the  acoustic  problem  may  need  more  accuracy  at  higher  angles,  for  the 
small  propagation  angles  important  in  the  electromagnetic  case,  Dockery  [14]  found  no  significant 
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impriyvement  between  the  SPE  propagator  and  the  WARE  propagator.  Therefore,  the  remainder  of  this 
section  will  be  dealing  with  the  approximation  given  by  Oj  (x)- 


0i  (x)  can  now  be  used  in  equation  13; 
ikr,Ax(A  +  B) 

u(x  +  Ax,  z)  =  e  u(x,z) 


(17) 


where 


A  =  (l/2)a;  B  =  (l/2)b. 

The  second-order  partial  derivative  in  the  exponent  can  easily  be  handled  by  using  Fourier  transforms. 
The  Fourier  transform  is  defined  as 

U(x,p)  =  .?'(u(x,z)l  =  f*  u(x,z) e'''P^dz  (18a) 

J  oo 

u(x,z)  =  .^■‘[U(x,p)]  =  _L  f*  U(x.p)e'P^  dz  (18b) 

2rr 


where  the  transform  variable  p  is  sin  0,  0  being  the  propagation  angle  above  the  horizontal.  Using 
the  transform  property  that 


~  (-ip)"  U(x,p) 


(19) 


az" 


and  noting  that  A  and  B  commute  (A  4-  B  =  B  -l-  A) 


ik  AxB 

u(x  +  Ax,z)  =  e 


ik„  Ax  A 


u(x,z) 


(20) 


equation  17  can  now  be  solved  to  give 


ik^:^(m--l) 

u(x  +  Ax,z)  =  e  "  .9' 


Ax  p- 


2  k, 


”  .?'|u(x.z)l 


(21) 


Equation  21  is  the  formula  that  must  be  computed  at  each  range  step.  Ax,  based  on  the  field  at  the 
previous  range. 

When  the  refractive  index,  m,  is  range  dependent,  the  quantities  A  and  B  no  longer  commute,  in 
which  case  the  steps  leading  to  the  right-hand  side  of  equation  13  become  more  complicated. 
Approximations  and  assumptions  have  to  be  made  in  which  one  finally  obtains  the  same  equation  for  the 
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field  (equation  21),  but  with  a  somewhat  complicated  error  term.  This  error  term  depends  on  the  range 
step  size.  Ax,  the  frequency,  and  the  refractive  index  gradients.  Therefore,  one  can  make  the  errors 
associated  with  a  range-dependent  m  in  equation  21  small  by  taking  a  sufficiently  small  step  size  and 
assuring  that  m(x,z)  varies  slowly  with  range. 

2.3  PROPAGATION  FACTOR  AND  PATH  LOSS 

The  final  step  is  to  determine  the  pattern  propagation  factor,  F,  once  the  field  u(x,z)  has  been 
computed.  The  pattern  propagation  factor  is  defined  as  the  ratio  of  the  magnitude  of  the  field  at  a  point 
in  space  lEl,  to  the  magnitude  of  the  field  at  the  same  point  under  free-space  conditions,  IEqI; 


(22) 


The  field  attributable  to  a  horizontally  polarized  source  in  free  space,  located  on  the  vertical  axis,  is 
proportional  to  (in  spherical  coordinates) 


ikoR 


R 


(23) 


R  =  (r^  •'o  " 

where  r^  is  the  distance  to  the  source  from  the  point  of  reference  and  0  is  the  angle  with  respect  to  the 
vertical  between  the  source  and  field  point.  Recalling  equation  1,  the  propagation  factor  is  now  written 


^<)>  _  lAlr(r,0)lR^ 
(r  sin0)^^^ 


(24) 


where  the  normalizing  constant  is  absorbed  in  li|i  (r,0)l.  Changing  from  spherical  to  Cartesian  coordinates, 
R  becomes 


R  =  [x2  +  (z-z/|'^2  ^  ^ 


and  the  denominator  becomes 


1  + 


(z  -  z„) 


1/2 


X  for  -  «:! 


(25) 


r  sin0  =  a  sin0  a  a0  a  x  for  x/a  «\ 


(26) 


to  obtain  for  F 

F  =  li|i(x,z)l/x'  =  lu(x,z)l  f/x  . 


(27) 
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u(x,z)  is  a  complex  function,  possessing  a  real  and  imaginary  component 

u(x,z)  =  +  i  Uy  (28) 

In  dB,  F  is  then 

10  log  F"  =  iO  log  [(u^  +  Uy)  X  |.  (29) 

Path  loss,  in  dB,  is  defined  as 

loss  =  -10  log  ^  .  (30) 

(4ji  x)- 

Using  equation  29,  path  loss  is  calculated  as 

9  9  9 

X-  (u  +  u') 

loss  =  -10  log  - - L  .  (31) 

(4ii)-x 

2.4  NUMERICAL  IMPLEMENTATION 

The  Fourier  transform  in  equation  21  is  implemented  by  using  an  FFT  algorithm.  The  transform 
size  is  kept  at  a  fixed  number.  The  maximum  propagating  angle  above  the  horizontal  and  the  maximum 
height  is  made  to  depend  on  the  frequency,  thereby  limiting  the  program  to  low-altitude,  small-angle 
propagation  effects  at  high  frequencies.  The  sample  size  is  based  on  Nyquist's  criterion 

^max  *  Pmax  =  (32) 

where  z^^^  is  the  maximum  height  in  the  calculation  domain  and  sinOj^^^,  0j^  being  the 

full  angular  width  above  and  below  the  horizontal  for  a  full  complex  transform.  Due  to  the  nature  of  the 
DFT,  the  field  is  abruptly  truncated  at  2^^^’  ‘‘'^suiting  in  reflections  from  the  nonphysical  upper  boundary. 
Filters  are  used  to  attenuate  the  field  smoothly  at  large  heights  and  large  propagation  angles,  to  keep 
reflections  from  entering  the  physical  portion. 

The  field  at  any  point  above  the  Earth  can  be  found  by  using  image  theory  and  the  appropriate 
boundary  conditions.  For  a  perfectly  conducting  surface  and  horizontal  polarization,  the  field  must  vanish 
at  the  surface:  u(x,0)  =  0.  The  next  .step  is  to  find  an  initial  field  lu(0,z)]  for  which  equation  21  can 
then  be  used  to  "march"  the  field  out  in  range.  One  begins  by  noting  that  the  field  at  range  zero  is 
essentially  the  antenna  aperture  distribution,  and  that  the  far-field  antenna  pattern  and  its  aperture 
distribution  are  a  Fourier  transform  pair; 

u(o,z)  U(o,p)  =  f(p).  (33) 

The  antenna  pattern  f(p)  is  assumed  to  be  symmetric  about  p  =  0  (0  elevation  angle).  Introducing  a 
source  height  and  elevation  angle  is  easily  done  by  using  the  Fourier  transform  shift  properties 


II 


u(z  -  Zj,) 


.9^  -'PZs 
— ►  e 


t(p). 


(34a) 


u(0,z) 


Po)' 


(34b) 


where  z^  is  the  source  height  and  is  sin  0gj,  being  the  elevation  angle. 

Zero  elevation  angle  was  used  in  all  results  presented  in  this  report.  A  uniform  aperture 
distribution  (essentially  an  omnidirectional  antenna  [f(p)  =  1  ])  was  used  for  the  antenna  pattern,  making 
the  initial  field  in  p-space 


-ipz«  ipz^. 

U(0,p)  =  f(p)e  -  f(-p)  e  ■  =  -2iC„  sin(pZs)  (35) 

where  Cj,  is  the  normalizing  constant  and  the  two  terms  on  the  right-hand  side  of  U(0,p)  represent  the  real 
and  image  fields,  respectively.  Note  that  U(0,p)  is  an  odd  function,  resulting  from  a  horizontally 
polarized  source  with  the  boundary  condition  for  a  perfectly  conducting  surface.  The  initial  field  for  a 
vertically  polarized  source  and  its  boundary  condition  (Bu/BzI^^q  =  0)  leads  to  an  even  function. 

A  troposcatter  model  is  included  in  PEPC  and  was  taken  from  Yeh  [15]  and  Rice,  et  al.  [16].  For 
a  full  description  of  the  model,  refer  to  NOSC  TD  1342  [17],  pp  1 16  120. 

To  keep  the  error  term  (associated  with  a  nonconstant  m  in  equation  21)  small,  we  must  keep  Ax 
small.  Ax  is  made  to  depend  on  the  wavelength,  thereby  forcing  the  program  to  take  small  range  steps 
at  large  frequencies  where  it  is  needed.  For  a  range-dependent  environment,  a  new  profile  is  obtained 
at  each  range  step  from  interpolation  of  user-specified  profiles.  The  method  of  interpolation  will  be 
discussed  in  section  3.2. 
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3.0  VALIDATION  OF  PE  MODEL 


3.1  COMPARISON  WITH  RANGE-INDEPENDENT  ENVIRONMENTS 

Some  comparisons  were  made  against  a  waveguide  program  called  MLAYER,  using  several  range- 
independent  environments.  MLAYER  was  developed  at  NOSC  by  Baumgartner  [18]  and  Pappert^ 
Figures  6  and  7  show  excellent  agreement  between  PEPC  and  MLAYER  for  a  300-meter  surface-based 
duct  at  L  band,  and  a  20-meter  evaporation  duct  at  X  band,  respectively.  Field  strength  in  dB  relative 
to  free  space  is  plotted  versus  height.  A  standard  atmosphere  at  S  band  was  also  used  for  comparison, 
but  results  between  the  two  models  were  indistinguishable,  so  that  case  is  not  shown  here. 

A  radiometeorological  investigation  was  conducted  jointly  by  the  Departments  of  Scientific  and 
Industrial  Research  of  the  L'nited  Kingdom  and  New  Zealand  under  the  name  of  the  Canterbury  Project 
in  Canterbury  Province  in  the  South  Island  of  New  Zealand  [10).  Transmitters  were  placed  in  aircraft 
that  flew  a  sawtooth  pattern  from  over  100  km  offshore  inward  toward  the  receiving  ground  stations 
located  at  Wakanui  Beach.  Data  were  taken  from  the  Canterbury  Project  for  the  afterncxtn  of 
5  August  1947.  This  case  was  chosen  becau.se  the  environment  varied  little  with  range.  A  strong 
surface  duct  was  present  from  20  km  offshore  out  to  100  km,  as  shown  by  figure  8.  Slant  paths  were 
used  for  the  model  to  represent  different  "legs"  of  the  aircraft  flight  paths.  Figures  9  and  10  show  the 
X-band  and  S-band  cases,  respectively,  against  the  observed  radio  data.  Both  figures  display  field 
strength  in  dB  relative  to  free  space  versus  height  for  one  "leg"  from  65  to  75  nautical  miles,  in  which 
the  aircraft  descended  from  1000  ft  (at  65  nmi)  to  near  the  ocean  surface  (at  75  nmi).  The  model 
compares  well  with  the  observed  data. 

3.2  COMPARISON  WITH  RANGE-DEPENDENT  ENVIRONMENTS 

Three  profiles,  along  with  their  corresponding  ranges,  were  specified  as  inputs  to  PEPC  to  simulate 
an  environment  that  varied  from  a  standard  atmosphere  at  the  antenna,  to  a  3(X)-meter  surface-based  duct 
at  1(X)  km,  and  back  to  a  standard  atmo.sphere  at  200  km.  Figure  1 1  .shows  a  lateral  heterogeneous 
raytrace  for  this  environment,  using  a  program  developed  at  NOSC  [2].  The  dotted  portion  of  the  plot 
shows  how  the  trapping  layer  varies  with  range  in  height  and  thickness.  The  PE  model  requires  a  profile 
at  each  range  step.  The  three  4-point  profiles  and  their  ranges  were  specified  as  inputs  and,  internal  to 
the  program,  linear  interpolation  was  used  to  obtain  the  height  and  M-unit  value  for  each  point  at  each 
range  .step  based  on  the  two  nearest  specifled  profiles.  This  method  of  interpolation  provides  for  logical 
and  smooth  transitioning  of  the  various  features  between  specified  profiles.  Figure  12  shows  the  PEPC 
coverage  diagram  for  this  environment.  Comparing  figures  1 1  and  12,  one  can  see  that  the  field  contours 
from  PEPC  follow  the  ray  paths  in  figure  1 1. 


*  Pappert,  R.  A.  1984.  "Field  Strength  and  Path-loss  in  a  Multilayer  Tropospheric  Waveguide  Environment," 
NOSC  TN  i.166  (October).  Naval  Ocean  Systems  Center.  San  Diego,  California.  NOSC  Technical  Notes 
are  working  documents  and  do  not  represent  an  official  policy  statement  or  the  Naval  Ocean  Systems  Center. 
For  further  information,  contact  the  author. 
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3.2.1  Guudulupe  Island  Measurements 

Radio  and  meteorological  data  were  measured  concurrently  in  over-water  paths  from  Guadalupe 
Island  to  San  Diego  during  1947  and  1948  [9],  Figures  13  through  15  show  the  data  records  for  3  days 
of  measurements.  In  this  experiment,  receivers  were  located  on  the  ground  and  the  transmitting  and 
meteorological  equipment  was  located  in  the  aircraft.  The  flight  paths  are  shown  in  the  upper  portion  of 
the  figures.  The  aircraft  took  radio  measurements  on  both  ascending  and  descending  paths,  but  only  made 
meteorological  measurements  on  the  ascending  paths.  An  attempt  was  made  to  represent  the  ’’slant” 
profiles  vertically  at  ranges  corresponding  to  the  optimum  coupling  height.  (The  optimum  coupling  height 
is  the  height  at  which  the  base  of  the  inversion  layer  is  measured,  and  at  which  placement  of  an  antenna 
will  yield  maximum  angular  trapping.) 

The  measurements  of  12  March  1948  show  a  slowly  increasing  trapping  layer  from  40  to  200  nmi. 
Five  linear  segments  were  taken  of  each  profile.  In  figure  16,  the  five  profiles  are  displayed  in  M-units 
at  each  range  corresponding  to  the  optimum  coupling  height.  The  change  in  thickness  and  height  of  the 
most  important  features  of  the  profiles  are  shown.  Compari.sons  between  PEPC  with  observed  data  are 
shown  in  figures  17  through  25  for  the  170-,  520-,  and  3300-MHz  cases.  Again,  slant  paths  were  used 
for  the  model  to  simulate  the  flight  paths  of  the  aircraft.  All  height-gain  plots  are  displayed  in  dB  relative 
to  free  space  versus  height.  Standard  atmosphere  is  shown  as  a  reference. 

At  170  and  520  MHz,  there  is  excellent  agreement  at  all  three  slant  paths  between  PEPC  and  the 
observed  data.  For  the  S-band  case,  figures  23  and  24  agree  well,  while  at  greater  ranges  (figure  25) 
PEPC  begins  to  deviate  at  low  heights  from  the  measured  data.  This  discrepancy  may  be  attributable  to 
the  undetected  presence  of  an  evaporation  duct,  which  would  affect  higher-frequency  signals.  Accounting 
for  an  evaporation  duct  would  give  higher  signal  levels  at  the  lower  receiver  heights.  Figures  26  through 
28  are  PEPC  coverage  diagrams  for  the  170-,  520-,  and  3300-MHz  cases,  respectively.  In  the  S-band 
case  (figure  28).  it  is  more  apparent  how  the  environment  changes.  Comparing  figures  16  and  28,  it  is 
easy  to  see  how  the  field  follows  the  trapping  layer  as  it  varies  with  range. 

To  a.s.sume  a  homogeneous  environment  will  not  give  a  clear  picture  of  what  is  really  happening. 
Figure  29  shows  a  coverage  diagram  for  a  homogenous  environment  at  S  band,  using  the  profile  at  39  nmi 
in  figure  16.  While  the  coverage  diagram  may  match  that  of  figure  28  at  close  ranges,  it  differs  greatly 
at  large  ranges.  One  may  try  to  ’’average"  the  profiles  and  assume  that,  at  all  ranges,  the  field  is  "close 
to"  the  field  under  actual  range-dependent  conditions,  but  using  the  "average"  profile  is  also  unsuitable. 
The  average  profile  was  taken  from  the  data  record  in  figure  13  for  two  height-gain  plots  for  170  and 
3300  MHz  (figures  30  and  31,  respectively).  The  difference  between  the  predicted  and  the  measured 
fields  are  substantial. 

Height-gain  plots  are  shown  in  figures  32  through  43  tor  the  environment  on  8  April  1948.  Again, 
there  is  good  agreement  between  PEPC  and  the  measured  data  at  low  frequencies.  At  S  band,  there  is 
a  big  discrepancy  for  large  ranges  and  high  receiver  elevations.  One  can  only  speculate  that  the  presence 
of  an  elevated  trapping  layer  went  undetected,  yielding  low  signal  levels  at  large  ranges.  The 
corresponding  coverage  diagrams  for  the  three  frequencies  are  shown  in  figures  44  through  46. 

Figures  47  through  55  are  the  height-gain  plots  for  the  environment  of  13  November  1947.  There 
is  excellent  agreement  at  all  frequencies  and  ranges.  The  coverage  diagrams  are  shown  in  figures  56 
through  58. 
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3.2.2  Canterbury  Project  Measurements 

Figures  59  and  60  show  two  days  of  refractivity  measurements  from  the  Canterbury  project.  The 
aircraft  flight  paths  are  shown  on  the  right-hand  side.  The  frequencies  used  in  this  experiment  were 
predominantly  S  and  X  band,  although  some  measurements  were  also  made  at  60  MHz.  Two  antenna 
heights  were  used  on  the  receiving  ground  station,  at  26.5  and  86  ft.  In  the  following  plots,  SI  channel 
will  refer  to  S  band  with  antenna  height  of  26.5  ft,  and  S2  channel  will  refer  to  S  band  with  antenna 
height  of  86  ft.  Similarly,  XI  and  X2  channels  will  refer  to  X  band  with  antenna  heights  of  26.5  and 
86  ft,  respectively. 

Figures  61  through  63  show  height-gain  plots  for  the  environment  measured  on  19  June  1947.  The 
slant  path  used  corresponded  to  flight  leg  M-N.  There  is  excellent  agreement  in  all  three  cases.  The 
corresponding  coverage  diagrams  are  shown  in  figures  64  through  66. 

For  the  environment  of  1 1  July  1947,  two  flight  legs  were  chosen.  Figures  67  and  68  .show  leg 
G-H  for  the  XI  and  X2  channels,  respectively,  and  figures  69  and  70  show  the  same  channels  for  leg 
H-I.  Good  agreement  is  shown  for  leg  G-H,  with  somewhat  good  agreement  for  leg  H-'I.  The 
corresponding  coverage  diagrams  are  shown  in  figures  71  and  72. 


360  n  Surface-Based  Duct 


Frequency  (MHz);  1300.  Receiver  rgs  (kn  ):  100,0 

Pntenna  height  (n  );  30.5 

Mind  speed  (kn)  ;  0. 

Figure  6.  Height-gain  plot  comparison  of  PEPC  and  M LAYER  for  300-m 
surface-ba.sed  duct. 
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HEIGHT  (m) 


1000 


0 


20  n  Evaporation  Duct 


Frequency  (MHz):  9600.  Receiver  rgs  (kn  ):  100.0 

Rntenna  height  (n  ):  30. S 

Hind  speed  (kn)  :  0. 


Figure  7.  Height-gain  plot  comparison  of  PEPC  and  MLAYER  for  20-m 
evaporation  duct. 


BEACH 


310  330 


20  km 


310  330 


ASSUMED  PROFILE 


M  CURVES 


Figure  8.  Measured  M-unit-versus-height  profiles  from  Canterbury 
Project  environmental  measurements  for  5  August  1947. 


height  It!,  height  I'ti 


Canterbury  Project  S  August  19*17  -  XI  Channel 


Frequency  (HHz)  :  9875.  Slant  path  fron  -  to 

Antenna  height  (ft):  26.5  ranges  (nni):  65.8-  75.8 

Mind  speed  (kn)  :  0.  heights(ft)  :  1888.8-  8.0 


Figure  9.  X-band  comparison  between  PEPC  and  Canterbury  Project 
measurements  for  5  August  1947. 


Canterbury  Project  5  August  19*17  -  SZ  Channel 


Frequency  (MHz);  32*18.  Slant  path  fron  -  to 

Antenna  height  (ft):  86.8  ranges  (nni):  65.0-  75.0 

Hind  speed  (kn)  :  0.  heights(ft)  1000.8-  0.0 


Figure  10.  S-band  comparison  lictween  PEPC  and  Canterbury  Project 
mea.surements  for  5  Augu.st  1947. 
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Figure  13.  Environment  and  radio  data  record  from  Guadalupe  Island  mea.surements  for  12  March  1948. 
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Figure  14.  Environment  and  radio  data  record  from  Guadalupe  Island  measurements  for  8  April  1948. 


DISTANCE  IN  NAUTICAL  MILES  FROM  SAN  DIEGO 
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Environment  and  radio  data  record  from  Guadalupe  Island  measurements  for  13  November  1947 


HEIGHT  (ft) 


GUADALUPE  ISLAND  PROFILES 
12  MARCH  1948 


3500 

2800 

2100 

1400 

700 

0 

330  450  330  450  330  450  330  450  330  450 

M-units 

39  nmi  85  5nmi  125  nml  160  nmi  193  nmi 

Figure  16.  Interpolated  M-unit-versus-height  profiles  for  fixed  ranges  from 
12  March  1948,  with  variation  in  trapping  layer  height  and  thickness  represented 
by  shaded  areas. 


Cuaddlupe  Island  -  12  March  1948 


Frequency  (MHz):  170.  Slant  path  fron  -  to 

Rntenna  height  (ft):  100.0  ranges  (nni):  69.0-  89.0 

Hind  speed  (kn)  :  0.  heights(ft)  :  3500.0-  0.0 


Figure  17.  Compari.sons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  12  March  1948,  along  slant  path  from  69  to  89  nmi. 
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HEIGHT  (ft) 


Cudddiupe  IsldiN)  -  12  ftdrch  1948 


Frequency  (MHz )  ;  170 , 

Hntenna  height  (ft);  100.0 
Hind  speed  (kn)  ;  0. 


Slant  path  fron  -  to 
ranges  (nni);  91.0-  101.0 
heights(ft):  0.0-3500.0 


Figure  18.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
tor  12  March  1948,  along  slant  path  from  91  to  101  nmi. 


Guadalupe  Island  -  12  March  194R 


Frequency  (MHz):  170. 
Rntenna  height  (ft):  108.0 
Hind  speed  (kn)  :  0. 


Slant  path  fron  -  to 
ranges  (nni):  IIS.B-  130.0 
heightstft):  3500.0-  0.0 


Figure  19.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  12  March  1948,  along  slant  path  from  1 15  to  130  nmi. 
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PROPAGATION  FACTOR  (dB) 


Frequency  (NHz):  S20.  Slant  path  fron  -  to 

ftntenna  height  (ft);  188.8  ranges  (nni);  69.8-  89.8 

Mind  speed  (kn)  ;  8.  heightslft) ;  3588.8-  8.8 

Figure  20.  Compari.sons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  12  March  1948,  along  slant  path  from  69  to  89  nmi. 


Guadalupe  island  -  12  March  1948 


PROPAGATION  FACTOR  (dB) 

Frequency  (MHz);  528.  Slant  path  fron  -  to 

Rntenna  height  (ft);  188.8  ranges  (nni):  91.8-  181.8 

Hind  speed  (kn)  ;  8.  heights(ft)  ;  8. 8-  3588.8 


Figure  21.  Comparisons  between  PEPC  and  mea.sured  radio  data  at  520  MHz 
for  12  March  1948,  along  slant  path  from  91  to  101  nmi. 


HEIGHT  (ft)  HEIGHT  (ft) 


Cudddlupe  Isidnd  -  12  Mdrch  1948 


Frequency  (HHz)  :  520.  Slant  path  fron  -  to 

Antenna  height  (ft):  180.0  ranges  (nni):  115.0-  138.0 

Hind  speed  (kn)  .-  0.  heights(ft):  3588.8-  8. 8 

Figure  22.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  12  March  1948,  along  slant  path  from  1 15  to  130  nmi. 


Guadalupe  Island  -  12  March  1948 


Frequency  (KHz):  3300.  Slant  path  fron  -  to 

Antenna  height  (ft):  108.8  ranges  (nni):  80.0-  90.8 

Hind  speed  (kn)  :  8.  heights(ft) :  1688.0-  8.0 


Figure  23.  Compan.sons  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  12  March  1948,  along  slant  path  from  80  to  90  nmi. 
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HEIGHT  (ftl  HEIGHT  (ft) 


Cudddlupe  Isldnd  -  IZ  Mdrch  1948 


Frequency  (HHz):  3300.  Sldnt  path  fron  -  to 

Rntennd  height  (ft):  100.8  ranges  (nni);  91.0-  96.0 

Mind  speed  (kn)  ;  0.  heights(ft):  0.0-  1608.6 


Figure  24.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz 
tor  12  March  1948,  along  dit'tcrent  slant  paths  (91  to  96  nmi). 


Frequency  (MHz):  3300  Slant  path  fron  -  to 

Rntenna  height  (ft):  100  0  ranges  (nni):  124.0-  130.0 

Hind  speed  (kn):  0,  heights(ft)  :  1680.8-  0.0 


T  igure  2.S.  (’omparisons  between  PEPC  and  mea.sured  radio  data  at  3300  MHz 
ti'f  12  March  1948.  along  slant  path  from  124  to  130  nmi. 
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HEIGHT  (ft) 


3500 


Guadalupe  Island  -  12  March  1948 


2800 

2100 

1400 

700 
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PEPC  Mcr:  1.0 

Freq  MHz  170. 
ftntm  ft  100.0 
Mind  hn  0. 


0  30  60  90  120  150 


RANGE  (nmi) 

LOSS:0H|ll0|ill5  120Hil25Biil3eHl35  140nil45  150  999  dB 


Figure  26.  ('overage  diagram  tor  environment  of  12  March  1948  at  170  MHz. 
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HEIGHT  (ft) 


Guadalupe  Island  -  12  Harch  1948 


PEPC  Uer:  1.0 

Freq  MHz  528. 
RntHt  ft  100.0 
Mind  Kn  0. 


0  30  60  90 

_  RANGE  (nmi) _ 

LOSS:0Hil28H|l25  130Hil35nil44Hl45 


120  150 

15eHl55  160  999  dB 


Figure  27.  Coverage  diagram  ior  environment  of  12  March  1948  at  520  MHz. 
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HEIGHT  (ft) 


e  38  68  SB  128  158 

RANGE  (nmi) 

L0SS:8Hil38Hil35  148Hl459|l5eHl55  168flil65  178  999  dB 

Figure  28.  Coverage  diagram  for  environment  of  12  March  1948  at  3300  MHz. 
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HEIGHT  (ft) 


Guadalupe  Island  Honogeneous  -  12  March  1948 


0  38  60  90  120  150 

RANGE  (nmi) 

LOSS:0Hil30Hil35  140Hil45Bil5flil55  160Mil65  170  999  dB 

Figure  2^  Coverage  diagram  using  homogeneous  profile  measured  at  39  nmi 
taken  from  the  12  March  1948  data  record. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Cudddlupe  island  Ruerage  Profile  -  12  March  1948 


Frequencg  (MHz):  170.  Slant  path  fron  -  to 

Rntenna  height  (ft):  100.8  ranges  (nni):  91.0-  101.0 

Hind  speed  (kn):  0.  heights(ft)  :  B.B-  3506. B 


Figure  30.  Comparisons  between  PEPC  and  measured  radio  data,  using  the 
average  profile  from  the  12  March  1948  data  record  for  170  MHz. 


Guadalupe  Island  fluerage  Profile  -  12  March  1948 


Frequency  (MHz):  3300.  Slant  path  fron  -  to 

Rntenna  height  (ft):  100.0  ranges  (noi);  91.0-  96.0 

Hind  speed  (kn)  :  0.  heights(ft)  ;  0.0-  1600.0 


Figure  31.  Comparisons  between  PEPC  and  measured  radio  data,  using  the 
average  profile  from  the  12  March  1948  data  record  for  3300  MHz. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Cudddlupe  Isidnd  -  8  April  1948 


Frequency  (HHz):  170.  Sldnt  pdth  fron  -  to 

Antenna  height  (ft):  168.0  ranges  (nni);  60.0-  72.0 

Mind  speed  (kn)  :  0,  heights(ft)  :  8.0-  4000.8 


Figure  32.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  8  April  1948,  along  slant  path  from  60  to  72  nmi. 


Guadalupe  Island  -  8  April  1948 


PROPAGATION  FACTOR  (dB) 

Frequency  (MHz):  170.  Slant  path  fron  -  to 

Antenna  height  (ft):  108.0  ranges  (nni):  81.0-  95.0 

Mind  Sr^-ad  (kn)  :  0.  heights(ft)  :  4008.0-  0.8 


Figure  33.  Compari.sons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  8  April  1948,  along  slant  path  from  81  to  95  nmi. 


Cuaddlupe  Island  -  8  April  1348 


Frequency  (MHz )  :  520 . 

Antenna  height  (ft):  100.0 
Hind  speed  (kn)  :  0. 


Slant  path  fron  -  to 
ranges  (nni) ;  68.0-  72.0 

heights(ft):  0.8-4000.0 


Figure  34.  Comparison.s  between  PEPC  and  measured  radio  data  at  520  MHz 
for  8  April  1948,  along  slant  path  from  60  to  72  nmi. 


Guadalupe  Island  -  8  April  1948 


Frequency  (MHz):  170.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.8  ranges  (nni):  118.0-  136.8 

Mind  speed  (kn):  0.  heights(ft):  4000.0-  0.0 


Figure  35.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  8  April  1948,  along  slant  path  from  1 18  to  136  nmi. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Cudddtupe  Istand  -  B  Rpril  1948 


Frequency  (MHz):  170.  Slant  path  from  -  to 

Hntenna  height  (ft):  180.8  ranges  (nni):  97.0-  107.8 

Hind  speed  (kn)  :  0.  heights(ft)  :  0.0-  4808.0 


Figure  36.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  8  April  1948,  along  slant  path  from  97  to  107  nmi. 


Guadalupe  Island  -  8  Rpril  1948 


Frequency  (MHz):  528.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.0  ranges  (nni);  81.0-  95.0 

Hind  speed  (kn)  :  0.  heights(ft):  4888.0-  8.0 


Figure  37.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  8  .April  1948,  along  slant  path  from  81  to  95  nmi. 
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Guadalupe  Island  -  8  April  1948 


Frequency  (MHz):  520.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.0  ranges  (nni):  97.0-  107.8 

Hind  speed  (kn)  :  0.  heights(ft):  0.0-  4006.0 


Figure  38.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  8  April  1948,  along  slant  path  from  97  to  107  nmi. 


Guadalupe  Island  -  8  April  1948 


Frequency  (MHz):  520.  Slant  path  fron  -  to 

Antenna  height  (ft):  108.0  ranges  (nni);  118.0-  136.8 

Hind  speed  (kn)  0.  heights(ft) :  4000.0-  0.0 


Figure  39.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  8  April  1948,  along  slant  path  from  118  to  136  nmi. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Guadalupe  Island  -  8  flpril  1948 


Frequency  (MHz):  3300.  Slant  path  fron  -  to 

Rntenna  height  (ft);  100.0  ranges  (nni):  60.0-  69.0 

Hind  speed  (kn)  ;  0.  heights(ft)  ;  0.0-  3000.0 


Figure  40.  Comparison.s  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  8  April  1948,  along  slant  path  from  60  to  69  nmi. 


Guadalupe  Island  -  8  April  1948 


Frequency  (MHz):  3300  Slant  path  froo  -  to 

Antenna  height  (ft):  100  0  ranges  (nni);  85.0-  95.0 

Hind  speed  (kn);  0.  heights(ft)  ;  3000.0-  0.0 


Figure  41.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  8  April  1948,  along  slant  path  from  85  to  95  nmi. 
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HEIGHT  (fl)  HEIGHT  (ft) 


Cudddlupe  Isidnd  -  8  April  1946 


Frequency  (MHz):  3300.  Slant  path  fron  -  to 

Antenna  height  (ft):  108.8  ranges  (nni):  97.0-  104.0 

Mind  speed  (kn)  :  0.  heights(ft>:  8.0-  3000.0 


Figure  42.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  8  April  1948,  along  slant  path  from  97  to  104  nmi. 


Guadalupe  Island  -  8  April  1948 


Frequency  (MHz):  3300.  Slant  path  fron  -  to 

Antenna  height  (ft);  108.0  ranges  (nni):  123.0-  136.0 

Hind  speed  (kn)  :  0.  heights(ft)  :  3000.0-  0.0 


Figure  43.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  8  April  1948,  along  slant  path  from  126  to  136  nmi. 
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Figure  45.  Coverage  diagram  for  environment  of  8  April  1948  at  520  MHz. 
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Guadalupe  Island  -  8  April  1948 
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PEK  Mer:  1.0 
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Figure  46.  Coverage  diagram  for  environment  of  8  April  1948  at  3300  MHz. 
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Cuddalupe  Island  -  13  Nowenber  1947 


Frequency  (NHz):  170.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.0  rai>ges  (nni):  41.5-  52.5 

Hind  speed  (kn)  :  0.  heights(ft)  :  0.0-  4080.0 


Figure  47.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  13  November  1947,  along  slant  path  from  41.5  to  52.5  nmi. 


Guadalupe  Island  -  13  Nowenber  1947 


Frequency  (MHz):  170.  Slant  path  frOB  -  to 

Antenna  height  (ft):  100  0  ranges  (nni):  56.0-  75.0 

Hind  speed  (kn)  :  0.  heightsift)  :  4000.0-  0.0 


Figure  48.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  13  November  1947,  along  .slant  path  from  56  to  75  nmi. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Cuaddlupe  Island  -  13  Nouenber  1947 


Frequency  (HH2):  170.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.0  ranges  (nni):  77.0-  90.0 

Hind  speed  (kn)  :  0.  heights(ft)  :  0.0-  4000.0 


Figure  49.  Comparisons  between  PEPC  and  measured  radio  data  at  170  MHz 
for  13  November  1947,  along  slant  path  from  77  to  90  nmi. 


Guadalupe  Island  -  13  Novenber  1947 


Frequency  (HH2)  ;  520.  Slant  path  fron  -  to 

Antenna  height  (ft);  100.0  ranges  (nni);  41.5-  52.5 

Hind  speed  (kn)  ;  0.  heights(ft):  0.0-  4006.0 


Figure  50.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  13  November  1947,  along  slant  path  from  41.5  to  52.5  nmi. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Guadalupe  Island  -  13  Nouenber  1947 


Frequency  (MHz):  520.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.8  ranges  (mti):  56.8-  75.0 

Hind  speed  (kn)  :  0.  heights(ft)  :  4008. 0-  0.0 


Figure  51.  Comparisons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  13  November  1947,  along  slant  path  from  56  to  75  nmi. 


Guadalupe  Island  -  13  Novenber  1947 


Frequency  (MHz):  520.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.0  ranges  (noi);  77.0-  90.0 

Hind  speed  (kn)  :  0.  heights(ft)  :  0.0-  4MO.0 


Figure  52.  Compari.sons  between  PEPC  and  measured  radio  data  at  520  MHz 
for  13  November  1947,  along  slant  path  from  77  to  90  nmi. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Cudddliipe  Isldnd  -  13  Nouenber  1947 


Frequency  (MHz):  330d . 
Hntennd  height  (ft):  100.0 
Hind  speed  (Itn)  :  0. 


Sldnt  pdth  fron  -  to 
ranges  (nni):  41.S-  50.0 

heights(ft):  0.0-3000.0 


Figure  53.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  13  November  1947,  along  slant  path  from  41.5  to  50  nmi. 


Cuddalupe  Island  -  13  Novenber  1947 


Frequency  (MHz):  3300.  Slant  path  fron  -  to 

Rntenna  height  (ft):  100.0  ranges  (nni);  61.0-  75.0 

Hind  speed  (kn)  ;  0.  heights(ft)  ;  3808.0-  8.8 


Figure  54.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  13  November  1947.  along  slant  path  from  61  to  75  nmi. 
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HEIGHT  (ft) 


Cuadd lupe  Island  -  13  Nouenber  1947 


PROPAGATION  FACTOR  (dB) 

Frequency  (HHz);  3300.  Slant  path  fron  -  to 

Antenna  height  (ft):  100.8  ranges  (nni):  77.0-  87.0 

Mind  speed  (kn)  ;  0.  heights(ft):  0.0-  3080.0 


Figure  55.  Comparisons  between  PEPC  and  measured  radio  data  at  3300  MHz 
for  13  November  1947,  along  slant  path  from  77  to  87  nmi. 


HEIGHT  (ft) 


Guadalupe  Island  -  13  Novenber  1947 


PEK  Uer:  1.0 

Freq  HHz  170 . 
PntHt  ft  100.0 
Nind  kn  0. 


0  30  60  90  1Z0  150 

RANGE  (nmi) 

LOSS:0Hill0Hill5  120Hl25^13eHl35  140^145  150  999  dB 


Figure  56.  Coverage  diagram  tor  environment  of  13  November  1947  at  170  MHz. 
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Figure  5  7.  Coverage  diagram  for  environment  of  13  November  1947  at  520  MHz. 
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HEIGHT  (ft) 


PETC  Ver:  1.0 

Freq  MHz  3380. 
RntHt  ft  100.0 
Mind  kn  0. 


0  30 


LOSS:0 


60  90  1Z8  150 

RANGE  (nmi) 

140Hil45Hjil5flil55  160Hil65  170 


999  dB 


Figure  58.  Coverage  diagram  for  environment  of  13  November  1947  at  3300  MHz. 
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ISOPLETHS  OF  MODIFIED  REFRACTIVE  INDEX  RADAR  /-jC  FLIGHT  FLAN  -ANSON  Rl 


HEIGHT  (ft)  HEIGHT  (ft) 


laee 


Canterbury  Project  19  June  1947  -  S2  Channel 


Frequency  (MHz):  3240.  Slant  path  fron  -  to 

Antenna  height  (ft):  86.0  ranges  (nni):  26.0-  38.0 

Mind  speed  (kn)  :  0.  heightsfft)  :  0.0-  1808.0 


Figure  61.  Comparison  between  PEPC  and  measured  radio  data  at  3240  MHz 
from  Canterbury  Project  data  record  for  19  June  1947. 


Canterbury  Project  19  June  1947  -  Ml  Channel 


PROPAGATION  FACTOR  (dB) 

Frequency  (MHz):  9675.  Slant  path  fron  -  to 

Antenna  height  (ft):  26.5  ranges  (nni):  26.0-  38.0 

Hind  speed  (kn)  :  0.  heights(ft):  0.0-  1008.0 


Figure  62.  Comparisons  between  PEPC  and  measured  radio  data  at  9875  MHz 
from  Canterbury  Project  data  record  for  19  June  1947,  with  antenna  height  at 
26.5  feet. 


HEIGHT  (ft) 


Canterbury  Project  19  June  1947  -  X2  Channel 


Frequency  (MHz);  9875.  Slant  path  fron  -  to 

Antenna  height  (ft);  86.8  ranges  (nni);  26.0-  38. 0 

Hind  speed  (kn);  8.  heights(ft):  0.0-  1000.0 


Fiyure  63.  Coriparisons  between  PEPC  and  measured  radio  data  at  9875  MHz 
troni  Canterbury  Project  data  record  for  19  June  1947,  with  antenna  height  at 
86  teet. 
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PEPC  Mer:  1.0 

Freq  NHz  3240 . 
flntHt  ft  86.0 
Nind  kn  0. 


155  160  999  dB 


environment  from  figure  61. 


HEIGHT  (ft) 


HEIGHT  (ft) 


5( 


height  (ft)  HEIGHT  (ft) 


1800 


Canterbury  Project  11  July  1947,  Leg  C->H 


Frequency  (MHz);  9875.  Slant  path  fron  -  to 

Hntenna  height  (ft);  26.5  ranges  (nni);  53.8-  65.0 

Hind  speed  (kn)  ;  0.  heights(ft):  1080. B-  0.0 


Figure  67.  Comparisons  between  PEPC  and  measured  radio  data  for  flight  path 
labeled  GH  from  Canterbury  Project  data  record  of  1 1  July  1947,  for  antenna 
height  at  26.5  feet. 


Canterbury  Project  11  July  1947,  Leg  C->H 


Frequency  (MHz);  9875.  Slant  path  fron  -  to 

Rntenna  height  (ft);  86.0  ranges  (nni);  53.0-  65.0 

Hind  speed  (kn) ;  0.  heights(ft):  1000.0-  0.0 

Figure  68.  Compari.s<;ns  between  PEPC  and  measured  radio  data  for  flight  path 
labeled  GH  from  Canterbury  Project  data  record  of  11  July  1947,  for  antenna 
height  at  86  feet. 
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HEIGHT  (ft)  HEIGHT  (ft) 


Cdnterbury  Project  11  July  1947,  Leg  H->1 


Frequency  (MH2)  :  9875.  Slant  path  fron  -  to 

Antenna  height  (ft);  26.5  ranges  (nni):  40.0-  53.0 

Mind  speed  (kn)  ■-  0.  heights(ft);  0.0-  1008.0 


Figure  69.  Comparisons  between  PEPC  and  measured  radio  data  for  flight  path 
labeled  HI  from  Canterbury  Project  data  record  of  1 1  July  1947,  for  antenna 
height  at  26.5  feet. 


Canterbury  Project  11  July  1947,  Leg  H->1 


Frequency  (MHz):  9075.  Slant  path  fron  -  to 

Antenna  height  (ft);  86.8  ranges  (nni);  40.0-  53.0 

Hind  speed  (kn)  ;  8.  heights(ft) ;  0.0-  1000.0 


Figure  70.  Compari.sons  between  PEPC  and  measured  radio  data  for  flight  path 
labeled  HI  from  Canterbury  Project  data  record  of  11  July  1947,  for  antenna 
height  at  86  feet. 
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Figure  71.  Coverage  di  igram  for  frequencies  and  environments  from  figures  67—70, 
antenna  height  at  26.5  feet. 
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1000 


Canterbury  Project  11  July  1947 >  X2  Channel 


Figure  72,  Coverage  diagram  tor  frequencies  and  environments  from  figures  67—70, 
antenna  height  at  86  feet. 
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4.0  CONCLUSIONS 


The  parabolic  equation  model,  PEPC,  has  been  shown  to  predict,  with  reasonably  good  accuracy, 
field  strengths  for  range-dependent  and  range-independent  environments.  While  improvements  must  be 
made  to  account  for  surface  roughness  and  other  environmental  effects,  for  the  simplest  case  used  here 
PEPC  agreed  quite  well  with  experimental  data. 

As  mentioned  before,  horizontal  homogeneity  occurs  in  the  troposphere  almost  86%  of  the  time. 
When  the  environment  is  inhomogeneous,  the  model  pre.sented  here  provides  a  means  by  which  the  field 
strengths  can  be  predicted  within  these  environments.  The  comparisons  given  in  this  report  between 
predicted  field  strengths  and  measured  data  within  measured  range-dependent  environments  show  the 
importance  of  including  horizontal  inhomogeneity  in  field  prediction  models.  Care  must  be  taken, 
however,  in  applying  meteorological  measurements  to  the  models.  A  detailed  knowledge  of  the 
refractivity  structure  involved  is  needed  to  make  the  best  possible  estimate  of  how  the  environment  is 
changing. 
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APPENDIX 


GUADALUPE  ISLAND  PROFILE  TABLES 


A-1 


Table  A-1.  M-unit  versus  height  profiles  and  corresponding 
ranges  from  Guadalupe  Island  data  record  tor  12  March  1948. 

Environment  for  12  March  1948 


Profile  #1  -  Range  at  0  nmi 
Height (ft)  M-units 

0.  337. 

540.  358.44 

803.407  324.736 

1217.17  334.888 

1231.  334.494 

3500.  447.400 


Profile  #2  -  Range  at  39  nmi 
Height(ft)  M-units 

0.  337. 

540.  358.44 

803.407  324.736 

1217.17  334.888 

1231.  334.494 

3500.  447.400 


Profile  #3  -  Range  at  85.5  nmi 
Height (ft)  M-units 

0.  337.169 

740.  365.34 

1080.63  343.702 

1490.61  356.039 

1652.63  351.889 

3500.00  430.919 


Profile  #4  -  Range  at  125.  nmi 
Height ( f t )  M-units 

0.  335.819 

1190.  382.84 

1574.55  357.962 

1889.38  371.111 

2096.94  371.011 

3500.00  430.957 


Profile  #5  -  Range  at  160.  nmi 
Height(ft)  M-units 

0.  333.676 

1145.00  378.220 

2140.00  370.715 

2584.38  398.294 

2923.08  404.423 

3500.00  429.802 


Profile  #6  -  Range  at  193.  nmi 
Height (ft)  M-units 

0.  331.056 

1420.00  382.920 

2387.50  383.544 

2718.44  410.029 

2892.62  399.773 

3500.00  430.849 


A-2 


Table  A-2.  M-unit  versus  height  profiles  and  corresponding 
ranges  from  Guadalupe  Island  data  record  for  8  April  1948. 


Environment  for  8  April  1948 


Profile  M  - 

Range  at  0  nmi 

Profile  #2  -  Range 

at  41.1  nmi 

Height{ ft) 

M-units 

Height { ft) 

M-units 

0.0 

340.007 

0.0 

340.007 

349.827 

352.781 

349.827 

352.781 

574.223 

309.753 

574.223 

309.753 

968.295 

317.703 

968.295 

317.703 

1207.64 

345.040 

1207.64 

345.040 

1488.62 

351.822 

1488.62 

351.822 

1503.41 

358.671 

1503.41 

358.671 

1809.71 

359.083 

1809.71 

359.083 

2194.20 

375.590 

2194.20 

375.590 

2293.07 

374.704 

2293.07 

374.704 

2563.72 

393.768 

2563.72 

393.768 

4000.00 

447.772 

4000.00 

447.772 

Profile  #3  - 

Range  at  56.9  nmi 

Profile  #4  -  Range 

at  93.2  nmi 

Height ( ft) 

M-units 

Height ( ft) 

M-units 

0.0 

341.728 

0.0 

341.806 

368.239 

352.042 

672.037 

368.492 

577.280 

308.902 

983.245 

320.817 

1005.35 

320.349 

1231.19 

343.194 

1260.99 

348.064 

1388.96 

347.223 

1475.73 

352.767 

1503.41 

358.671 

1503.41 

358.671 

1809.71 

359.083 

1809.71 

359.083 

2194.20 

375.590 

2194.20 

375.590 

2293.07 

374.704 

2293.07 

374.704 

2563.72 

393.768 

2563.72 

393.768 

4000.00 

447.772 

4000.00 

447.772 

Profile  #5  - 

Range  at  132.7  nmi 

Height ( ft) 

M-units 

0.0 

342.668 

733.564 

369.305 

1059.78 

323.572 

1251.12 

345.108 

1539.68 

361.967 

2030.24 

368.838 

2153.73 

373.425 

2182.62 

371.732 

2363.42 

384.754 

4000.00 

443.040 

A-3 


Table  A-3.  M-unit  versus  height  profiles  and  corresponding 
ranges  from  Guadalupe  Island  data  record  for  13  November  1947. 


Environment  for  13  November  1947 

Profile  #1  - 

Range  at  0  nmi 

Profile  #2  -  Range 

Height ( ft) 

M-units 

Height ( ft) 

0. 

317.9 

0. 

216.865 

314.73 

795.918 

1617.62 

362.87 

1116.23 

3084.17 

420.63 

1553.31 

4000. 

456.18 

3182.06 

4000.00 

Profile  #3  - 

Range  at  71.4  nmi 

Profile  #4  -  Range 

Height{ ft) 

M-units 

Height { ft) 

0.0 

331.331 

0.0 

393.097 

335.389 

518.528 

618.075 

351.704 

790.755 

820.575 

355.378 

926.370 

1008.75 

362.805 

1214.22 

1359.41 

348.897 

2017.49 

2155.46 

371.392 

2556.49 

3182.06 

414.038 

3684.10 

4000.00 

440.863 

4000.00 

at  37.5  nmi 
M-units 
322.762 
345.177 
340.451 
354.080 
414.038 
440.863 

at  126.9  nini 
M-units 
306.607 
321.198 
336.550 
335.692 
348.295 
373.394 
374.805 
410.625 
425.710 


Table  A-4.  M-unit  versus  height  profiles  and  corresponding 
ranges  from  Canterbury  Project  data  record  for  11  July  1947. 


Environment  for  11  July  1947  from  Canterbury  Project 


Profile  #1  - 

Range  at  0  nmi 

Profile  #2  -  Range 

at  10.8  nmi 

Height ( ft) 

M-units 

Height ( ft) 

M-units 

0. 

322.5 

0. 

335.5 

140. 

324. 

60. 

327. 

530. 

342.5 

200. 

330. 

1000. 

351. 

270. 

329. 

1000. 

350. 

Profile  #3  - 

Range  at  21.6  nmi 

Profile  #4  -  Range 

at  32.4  nmi 

Height ( ft) 

M-units 

Height ( ft) 

M-units 

0. 

337.5 

0. 

338. 

25. 

333.5 

40. 

334. 

105. 

336. 

165. 

339. 

200. 

327. 

285. 

327. 

1000. 

350. 

1000. 

350. 

Profile  #5  - 

Range  at  43.2  nmi 

Profile  #6  -  Range 

at  56  nmi 

Height ( ft) 

M-units 

Height! ft) 

M-units 

0. 

339. 

0. 

336. 

150. 

331. 

80. 

333. 

195. 

326. 

200. 

333. 

350. 

330. 

270. 

327. 

1000. 

351. 

1000. 

347.5 

A-5 


Table  A-5.  M-unit  versus  height  profiles  and  correspionding 
ranges  from  Guadalupe  Island  data  record  for  19  June  1947. 

Environment  for  19  June  1947 

Profile  #1  -  Range  at  0  nmi  Profile  #2  -  Range  at  10.8  nmi 


Height ( ft) 

M-units 

Height ( ft) 

M-units 

0. 

304. 

0. 

315. 

220. 

307. 

10. 

309. 

515. 

315.5 

130. 

314. 

1000. 

318.1 

195. 

306. 

405. 

323. 

500. 

317.5 

1000. 

335.8 

Profile  #3  -  Range  at  21.6  nmi  Profile  #4  -  Range  at  32.4  nmi 


Height (ft)  M- 

-units 

Height ( ft) 

M-units 

0. 

315. 

0. 

320. 

25. 

312. 

20. 

313. 

345. 

325. 

220. 

322. 

425. 

323.5 

330. 

322.5 

520. 

326. 

500. 

329.5 

600. 

325. 

730. 

336. 

1000. 

339.6 

1000. 

347. 

Profile  #5  -  Range 

at  43.2  nmi 

Height (ft)  M- 

■units 

0. 

307. 

375. 

320. 

480. 

322. 

650. 

330. 

780. 

328. 

1000. 

339. 
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